V 



Doc Ref. FP1 
Appl. No. 10/073,154 



WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




PCT 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification 6 
H04L 25/03, 27/26 



Al 



(11) International Publication Number: WO 99/52250 

(43) International Publication Date: 14 October 1999 (14.10.99) 



(21) International Application Number: PCT/US99/07422 

(22) International Filing Date: 5 April 1999 (05.04.99) 



(30) Priority Data: 

09/054,468 
60/087,336 
09/233,914 



3 April 1998 (03.04.98) US 
29 May 1998 (29.05.98) US 
21 January 1999 (21.01.99) US 



(63) Related by Continuation (CON) or Continuation-in-Part 
(CIP) to Earlier Application 

US 09/233,914 (CON) 

Filed on 2 1 January 1 999 (2 1 .0 1 .99) 



(71) Applicant (for all designated States except US): TELLABS OP- 

ERATIONS, INC. [US/US]; 4951 Indiana Avenue, MS-16, 
Lisle, IL 60532 (US). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): HARIKUMAR, Gopal 
[IN/IN]; 1705 Enchanted Forest, South Bend, IN 46637 
(US). MARCHOK, Daniel [US/US]; 1 4984 West Clear Lake 
Road, Buchanan, MI 49107 (US). RUDOFSKI, Kenneth, J. 
[US/US]; Unit 4, 737 West Aldine, Chicago, IL 60657 (US). 



(74) Agent: WISCHHUSEN, Carl, B.; Fish & Richardson P.C., 601 
Thirteenth Street, N.W., Washington, DC 20005 (US). 



(81) Designated States: AE, AL, AM, AT, AU, AZ, BA, BB, BG, 
BR, BY, CA, CH, CN, CU, CZ, DE, DK, EE, ES, FI, GB, 
GD, GE, GH, GM, HR, HU, ID, IL, IN, IS, JP, KE, KG, 
KP, KR, KZ, LC, LK, LR, LS, LT, LU, LV, MD, MG, MK, 
MN, MW, MX, NO, NZ, PL, PT, RO, RU, SD, SE, SG, SI, 
SK, SL, TJ, TM, TR, TT, UA, UG, US, UZ, VN, YU, ZA, 
ZW, ARIPO patent (GH, GM, KE, LS, MW, SD, SL, SZ, 
UG, ZW), Eurasian patent (AM, AZ, BY, KG, KZ, MD, 
RU, TJ, TM), European patent (AT, BE, CH, CY, DE, DK, 
ES, FI, FR, GB, GR, IE, IT, LU, MC, NL, PT, SE), OAPI 
patent (BF, BJ, CF, CG, CI, CM, GA, GN, GW, ML, MR, 
NE, SN, TD, TG). 



Published 

With international search report. 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. 



(54) Title: FILTER FOR IMPULSE RESPONSE SHORTENING, WITH ADDITION SPECTRAL CONSTRAINTS, FOR MULTICAR- 
RIER TRANSMISSION 



INPUT 
FROM _ 
RECEIVE 
FILTER 



AID 



yk 



100 



SOSF 



\ 



x k 



STRIP 
CP 



SCISF 
COPY 



\ 



yk 

r 30S 



DELAY 
BUOCK 



4 



205 



ADAPTATION 

\ 



yk 



r110 



'120 



-130 rf40 



SERIAL 

TO 
PARALLEL 



OFT 



FEQ 



OUTPUT 
; DECODER '"DATA 



ADD 
CP 



r21S 



IMS 
ALGORITHM 



V250 



PARALLEL 

TO 
SERIAL 



IDFT 



^■340 



ENCODER 
T 



(57) Abstract 

A channel in a multiple carrier communication system is equalized by computing a desired spectral response, shortening the impulse 
response of the channel so that a significant part of an energy of the impulse response is confined to a region that is shorter than a target 
length and filtering the signal based on the desired spectral response. A multiple carrier communication system may include a primary 
impulse shortening filter that receives an output signal of an analog to digital converter and accepts coefficients. A secondary impulse 
shortening filter may receive the output signal of the analog to digital converter, output an output signal, and pass coefficients to the primary 
impulse shortening filter. A reference signal generator may output a reference signal. A comparator may compare the output signal and the 
reference signal and output a resulting error signal. An adaptive processor may compute coefficients for the secondary impulse shortening 
filter based on the error signal. 
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FILTER FOR IMPULSE RESPONSE SHORTENING, WITH ADDITION SPECTRAL CONSTRAINTS, 
FOR MULT I CARRIER TRANSMISSION iwiiuuma. 



Background 

5 The invention relates to time-domain equalization in a discrete multi-tone (DMT) 

receiver. 

Conventional single carrier modulation techniques translate data bits for transmission 
through a communication channel by varying the amplitude and/or phase of a single sinusoidal 
carrier. By contrast, DMT, which is also referred to as Orthogonal Frequency Division 
1 0 Multiplexing (OFDM) or Multicarrier Modulation (MCM), employs a large number of sinusoidal 
subcarriers, e.g., 128 or 256 subcarriers. The available bandwidth of the communication channel 
is divided into subchannels and each subchannel communicates a part of the data. A DMT 
system may employ quadrature amplitude modulation (QAM) for each of the subcarriers. 

OFDM-based systems transmit blocks of information bits. The time required to transmit 
15 one such block is called the symbol period. The time domain waveform that corresponds to one 
such block of bits is called a symbol. 

Intersymbol interference (ISI) arises from the characteristics of practical communication 
channels and limits the rate at which information can be transmitted through them. Specifically, 
communication channels typically have an Effective Discrete-Time Impulse Response (EDER) 

2 0 that is greater than one sample time in length, which causes ISI. ISI is a well-known 

phenomenon in single-carrier communication systems and there are many techniques for 
reducing it. The process of such ISI reduction is called equalization. ISI is discussed, for 
example, in Proakis, Digital Communications, McGraw Hill, 2nd Edition, 1989. 

Equalization in OFDM-based systems is achieved by a two stage process. First, at the 
25 transmitter, a Cyclic Prefix (CP) is employed by affixing an end-portion of each symbol to the 
beginning of the symbol. A cyclic prefix that is greater than the EDIR of the channel prevents 
one symbol from interfering with another. Furthermore, it also facilitates a simple method of 
neutralizing the time domain spread of each symbol forced by the channel. This is achieved 
through a simple frequency domain process in the receiver which requires one multiplication 

3 0 operation for each subcarrier used. The use of a Cyclic Prefix to reduce ISI is discussed, for 

example, in: Cimini, "Analysis and Simulation of a Digital Mobile Channel using Orthogonal 
Frequency Division Multiplexing," IEEE Transactions on communications, pp 665-675, July 
1985; Chow, "A Discrete Multi-Tone Transceiver System for HDSL applications," IEEE Journal 
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on Selected Areas of Communications, 9(6):895-908, August 1991; and "DMT Group VDSL 
PMD Draft Standard Proposal," Technical Report, T1E1.4/96-329R2, ANSI 1997. 

Another problem arising in conventional DMT systems is noise bleeding, which occurs 
when noise from one frequency band interferes with a signal whose subcarrier is in another 
5 frequency band. Noise bleeding is caused, in general, by a discrete Fourier transform (DFT) 
operation at the receiver. Noise bleeding is discussed in, for example, Worthen et. aL, 
"Simulation of VDSL Test Loops," Technical Report T1E1. 4/97-288, ANSI 1997. 

In a perfectly synchronized DMT system, a signal in one frequency band does not 
interfere with a signal whose subcarrier is in another frequency band. However, noise from one 
1 0 band may interfere with other less noisy bands and render them unusable. Techniques for 

dealing with noise-bleeding include wavelet-based solutions. However, wavelet-based solutions 
are, in general, computationally intensive. 

Other references dealing with time domain equalization include: Chow, J.S. and Cioffi, 
J.M., "A Cost-effective Maximum Likelihood Receiver for Multicarrier Systems", Proceeding s 
1 5 of the ICC . 1 992; Melsa, Peter J.W., Younce, Richard C, and Rohrs, Charles E., "Optimal 

Impulse Response Shortening", Proceedings of the thirty-third Annual Allerton Conference on 
Communication. Control and Computing. 1995, pp. 431-438; Harikumar, Gopal and Marchok, 
Daniel, "Shortening the Channel Impulse Response of VDSL Loops for Multicarrier 
Applications", Technical report T1E1. 4/97-289, ANSI, 1997. 

20 

Summary 

The invention provides a spectrally constrained impulse shortening filter (SCISF) for use, 
for example, in DMT systems. The SCISF serves two primary functions. 

First, the SCISF reduces intersymbol interference (ISI) by reducing the length of the 

2 5 effective discrete-time impulse response (EDIR) of the communication channel. Conventional 

impulse shortening filters may have deep nulls in their frequency response. By contrast, the 
SCISF has a filter characteristic that is essentially free from undesired nulls that may attenuate or 
completely eliminate certain subcarriers. 

Second, the SCISF reduces noise bleeding between subchannels by attenuating noisy 

3 0 channels in a manner that does not reduce the signal to noise ratio (SNR) in these channels, but 

reduces the noise power that may appear in the sidelobes of adjacent subchannels. The SCISF 
accomplishes these functions by applying a frequency constraint to the signal based on a desired 
spectral response. 
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In one general aspect, the invention features equalizing a channel in a multiple carrier 
communication system. The system includes a spectrally constrained impulse shortening filter. 
Received noise power spectral density is measured. A desired spectral response is computed 
based on the measured noise power. A frequency response for the spectrally constrained impulse 
5 shortening filter is selected based on the desired spectral response. The communication signal is 
filtered with the spectrally constrained impulse shortening filter. 

Embodiments may include one or more of the following features. The noise power 
spectral density may be measured at the output of the discrete Fourier transform. The spectrally 
constrained impulse shortening filter may be a time domain digital filter. 
10 In another aspect, the invention features equalizing a channel in a multiple carrier 

communication system. The channel has an impulse response and is configured to receive a 
signal having a cyclic prefix. A desired spectral response is computed. The impulse response of 
the channel is shortened so that a significant part of the energy of the impulse response is 
confined to a region that is shorter than a target length, and the signal is filtered based on the 
1 5 desired spectral response. 

The target length may be a length of the cyclic prefix. The desired spectral response may 
be computed from measured noise power density, which may be measured at the output of a 
discrete Fourier transform. For example, the desired spectral response may be the inverse of the 
measured noise power spectral density. The filtering step may be performed with a filter having 
2 0 a frequency response selected to match the desired spectral response. The shortening of the 
impulse response and/or the filtering may be performed by a time domain digital filter. 

In another general aspect, the invention features selecting an impulse response for a 
spectrally constrained impulse shortening filter in a multiple carrier communication system. 
Received noise power spectral density is measured. A cost function is computed using the noise 

2 5 power. The cost function is dependent on the impulse response. The dimensionality of a space 

over which the cost function is defined is reduced and the cost function is minimized. The noise 
power spectral density may be measured at the output of a discrete Fourier transform. The cost 
function may be used to compute coefficients for the spectrally constrained impulse shortening 
filter. 

3 0 In another general aspect, the invention features a spectrally constrained impulse 

shortening filter for a multiple carrier communication system. The system includes a channel 
that has an impulse response. A digital filter structure receives the signal and apply a frequency 
characteristic to the signal. The frequency characteristic is determined by filter coefficients. 
Filter coefficients are selected to shorten the impulse response of the channel so that a significant 
3 5 part of the energy of the impulse response is confined to a region that is shorter than a target 
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length and to apply a frequency characteristic to the signal based on a desired spectral response. 
The selected filter coefficients are input to the taps of the filter. 

In another general aspect, the invention features a receiver for receiving a multiple carrier 
signal from a communication channel having an impulse response. An analog-to-digital 
5 converter receives the signal from the communication channel, A spectrally constrained impulse 
shortening filter receives the signal from the analog-to-digital converter and shortens the impulse 
response of the channel so that a significant part of an energy of the impulse response is confined 
to a region that is shorter than a target length. The filter also applies a frequency characteristic to 
the signal based on a desired spectral response. A discrete Fourier transform receives the output 

1 0 of the spectrally constrained impulse shortening filter and a decoder receives outputs of the 
discrete Fourier transform. 

In another aspect, the invention features a modem in which an encoder receives digital 
data and outputs a constellation of complex values. An inverse discrete Fourier transform 
receives the constellation from the encoder. A digital-to-analog converter is connected to the 

15 inverse discrete Fourier transform and outputs a signal to a communication channel. An analog- 
to-digital converter receives the signal from the communication channel. A spectrally 
constrained impulse shortening filter shortens an impulse response of the channel so that a 
significant part of an energy of the impulse response is confined to a region that is shorter than a 
target length. The signal is filtered based on a desired spectral response. A discrete Fourier 

2 0 transform connected to the filter and a decoder is connected to the discrete Fourier transform and 
outputs digital data. 

In another aspect, the invention features software for causing a processor in a 
communication system to perform the following operations: measure received noise power 
spectral density; and compute a desired spectral response based on the measured noise power. 

2 5 The software may also include instructions for causing a processor to compute filter coefficients 

based on the desired spectral response. 

In another aspect, the invention features software for causing a processor in a 
communication system to measure received noise power spectral density. A cost function is 
computed using the noise power. The cost function is dependent on an impulse response of a 

3 0 spectrally constrained impulse shortening filter. The dimensionality of a space over which the 

cost function is defined is reduced, and the cost function is minimized. 

The techniques described here are not limited to any particular hardware or software 
configuration. They may find applicability in any computing or processing environment that 
may be used for a communication system. The techniques may be implemented in hardware or 
3 5 software, or a combination of the two. Preferably, the techniques are implemented in computer 



BNSDOCID: <WO 9952250A1J_> 



WO 99/52250 PCT/US99/07422 

5 

programs executing on a digital signal processor that includes a processor and a storage medium 
readable by the processor (including volatile and non-volatile memory). 

The coefficients of the SCISF may be computed and changed through a training process 
that may occur upon start up of the communication system and periodically during its operation. 
5 Using the training process, the coefficients of the SCISF are computed independently of the 
length of the cyclic prefix, the symbol length, and the frequency domain equalization 
characteristics of the system. 

The SCISF is particularly effective in systems in which additive noise dominates 
intersymbol interference, or in which noise bleeding predominates. The SCISF reduces noise 

1 0 bleeding with a filter structure that is shorter than that obtained with other techniques. 

Consequently, the SCISF is less complex and its coefficients are easier to compute, which 
reduces system complexity and cost. The SCISF may be particularly well suited, for example, 
for very high-speed digital subscriber lines (VDSL) systems, which generally have low 
intersymbol interference and tend to suffer from noise bleeding. 

15 In addition, dynamic selection of a cyclic prefix (CP) that maximizes the data throughput 

for a communication channel having a particular noise profile is provided. Dynamic selection of 
the CP allows the communication system to adapt to changing noise conditions in the 
communication channel. 

In one aspect, generally, a primary impulse shortening filter is adapted in a multiple 

2 0 carrier communication system. A secondary impulse shortening filter is provided. An output 

signal of the secondary impulse shortening filter is compared to a reference signal to compute an 
error signal. Coefficients of the secondary impulse shortening filter are computed in an adaptive 
processor based on the error signal. Coefficients of the primary impulse shortening filter are 
replaced with coefficients of the secondary impulse shortening filter. 

2 5 Embodiments may include one or more of the following features. An output signal of the 

primary impulse shortening filter may be decoded to form output data. The output data may be 
encoded to form the reference signal. A discrete Fourier transform may be applied to the output 
signal of the primary impulse shortening filter prior to decoding the output signal. An inverse 
discrete Fourier transform may be applied to the encoded output data in forming the reference 

3 0 signal. 

A digital signal may be received from an output of an analog to digital converter. The 
digital signal may be input to the primary impulse shortening filter. The digital signal may be 
delayed. The delayed digital signal may be input to the secondary impulse shortening filter and 
the adaptive processor. 
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The encoded output data may be scaled with a set of scaling factors in forming the 
reference signal. The scaling factors may be determined by: measuring received noise power 
spectral density, computing a desired spectral response based on the measured noise power, and 
computing the scaling factors so that the coefficients computed in the adaptive processor provide 
5 the secondary impulse shortening filter with a spectral response that matches the desired spectral 
response. A discrete Fourier transform may be applied to the output signal of the primary 
impulse shortening filter prior to decoding the output signal. The noise power spectral density 
may be measured at an output of the discrete Fourier transform. An inverse discrete Fourier 
transform may be applied to the scaled, encoded output data. 

10 In another aspect, an impulse shortening filter is adapted in a multiple carrier 

communication system having a spectrally constrained impulse shortening filter. An output 
signal of the spectrally constrained impulse shortening filter is compared to a reference signal to 
compute an error signal. Coefficients of the spectrally constrained impulse shortening filter are 
computed in an adaptive processor based on the error signal. 

15 Embodiments may include one or more of the following features. The reference signal 

may be a predetermined signal stored in a memory in the communication system. A discrete 
Fourier transform may be applied to predetermined reference values to form transformed 
reference values. The transformed reference values may be scaled with a set of scaling factors to 
form scaled. An inverse discrete Fourier transform may be applied to the scaled values to form 

2 0 the reference signal. 

A data signal may be received from an output of an analog to digital converter. The data 
signal may be input to the spectrally constrained impulse shortening filter and the adaptive 
processor. 

In another aspect, a multiple carrier communication system includes a primary impulse 

2 5 shortening filter that receives an output signal of an analog to digital converter and accepts 

coefficients. A secondary impulse shortening filter receives the output signal of the analog to 
digital converter, outputs an output signal, and passes coefficients to the primary impulse 
shortening filter. A reference signal generator outputs a reference signal. A comparator 
compares the output signal and the reference signal and outputs a resulting error signal. An 

3 0 adaptive processor computes coefficients for the secondary impulse shortening filter based on the 

error signal. 

Embodiments may include one or more of the following features. A discrete Fourier 
transform may receive an output signal of the primary impulse shortening filter. A decoder may 
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receive the transformed output signal from the discrete Fourier transform. The reference signal 
generator may include an encoder that receives output data from the decoder. The reference 
signal generator may include a scaling filter that scales the output data from the encoder using a 
set of scaling factors. An inverse discrete Fourier transform may receive the scaled output signal 
5 from the scaling filer. 

In another aspect, a multiple carrier communication system may include a spectrally 
constrained impulse shortening filter that receives an output signal of an analog to digital 
converter and accepts coefficients. A reference signal generator outputs a reference signal. A 
comparator compares the output signal and the reference signal and outputs a resulting error 
1 0 signal. An adaptive processor computes coefficients for the spectrally constrained impulse 
shortening filter based on the error signal. 

Embodiments may include one or more of the following features. A memory may store 
the reference signal as a predetermined signal. A discrete Fourier transform may receive the 
reference signal from the memory. A scaling filter may scale the reference signal using a set of 
1 5 scaling factors. An inverse discrete Fourier transform may receive the scaled reference signal. 

Other features and advantages will be apparent from the following description, including 
the drawings and the claims. 

Brief Descripti on of the Drawings 
2 0 Fig. 1 is a block diagram of a discrete multi-tone communication system. 

Fig. 2 is a plot of effective discrete-time impulse response (EDIR) of a communication 
channel including transmit and receive filters. 

Fig. 3 is a plot of the shortened EDIR due to a spectrally constrained impulse shortening 
filter (SCISF). 

2 5 Fig. 4 is a block diagram of a SCISF. 

Fig. 5 is a plot of transmit signal power, signal power at SCISF input and noise power at 
SCISF input. 

Fig. 6 is a plot of signal and noise power at the output of the SCISF. 

Fig. 7 shows the filter response of a discrete Fourier transform for one frequency band. 

3 0 Fig. 8 is a plot of the desired spectral response of the SCISF, G/co), versus the actual 

frequency response, G(co). 

Fig. 9 is a plot of signal-to-noise ratio at the output of the SCISF, the output of the DFT 
with the SCISF and the output of the DFT without the SCISF. 
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Fig. 10 is a block diagram of a system for adapting a spectrally constrained impulse 
shortening filter in an initial training process. 

Figs. 1 1 and 12 are block diagrams of a system for adapting a spectrally constrained 
impulse shortening filter in an periodic or continuous training process. 
5 Fig. 13 is a block diagram of a generalized adaptation system that employs frequency 

scaling in the feedback loop. 

Fig. 14 is a block diagram of a system for adapting a spectrally constrained impulse 
shortening filter in an initial training process that includes frequency scaling in the feedback 
loop. 

10 Figs. 15 and 16 are block diagrams of a system for adapting a spectrally constrained 

impulse shortening filter in an periodic or continuous training process that includes frequency 
scaling in the feedback loop. 

Figs. 17-20 are plots of simulation results for a discrete multi-tone system. 

Figs. 21-24 are plots of simulation results for a discrete multi-tone system. 

15 

Description 

As shown in Fig. 1, a discrete multi-tone (DMT) communication system 10 has a 
transmitter 12 and a receiver 14. The transmitter 12 accepts an input data bit stream which 
passes through a constellation encoder 20. The encoder 20 divides the serial input bit stream into 
2 0 blocks of data. These blocks of data are further subdivided into smaller blocks corresponding to 
subchannels. Each of these smaller blocks are used to compute a complex value representing a 
constellation point. Each constellation point corresponds to a subsymbol. The subsymbols are 
then output by the encoder. Taken together, the subsymbols constitute a symbol. 

The subsymbols are supplied to an inverse discrete Fourier transform (IDFT) 30, which 

2 5 may be implemented, for example, in a digital signal processor. The IDFT 30 outputs N time 

samples of a symbol. The time samples are processed by a parallel to serial converter 40 to form 
a single stream of time samples. 

Following the parallel to serial converter 40, a prefix adder 50 adds a cyclic prefix to the 
beginning of each symbol to reduce intersymbol interference (ISI). Alternatively, the cyclic 

3 0 prefix may be added in the parallel to serial converter. After the cyclic prefix is added, the 

resulting signal passes through a digital-to-analog (D/A) converter 60 for transmission to the 
receiver 14 through a communication channel 70. An analog transmit filter 65 may be included 
following the D/A converter to band limit the transmitted signal. 
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At the receiver 14, the signal passes through an analog-to-digital (A/D) converter 80 and 
then through a spectrally constrained impulse shortening filter (SCISF) 90. A prefix stripper 100 
strips the cyclic prefixes from the resulting symbols and a serial to parallel converter 110 divides 
the stream of time samples into parallel signal paths that form the inputs to a discrete Fourier 
5 transform (DFT) 120. The DFT 120 converts the time samples into subsymbols. A frequency 
domain equalization filter 130 equalizes the subsymbols. A decoder 140 converts the 
subsymbols into a data bits and outputs the resulting data. An analog receive filter 75 may be 
included prior to the A/D converter in order to band limit the received signal. 

As discussed above, a cyclic prefix is added to each symbol prior to transmission through 

1 0 the communication channel to reduce the effects of ISI. The cyclic prefix is added by copying 
the last v time samples from the end of a symbol and placing them at the beginning of the 
symbol. To eliminate ISI, the length of the cyclic prefix, v, is chosen to be longer than the 
effective discrete-time impulse response (EDIR) of the channel. However, because the cyclic 
prefix constitutes redundant data, increasing the length of the cyclic prefix reduces the efficiency 

1 5 of the communication system. For example, in a system having N time samples per symbol and 
a cyclic prefix of v time samples, the efficiency of the system will be reduced by a factor of 
N/(N+v). Efficiency may be maximized either by minimizing v or maximizing N. However, 
increasing N increases the complexity, latency and computational requirements of the system and 
at some point becomes impractical. Accordingly, it is desirable to minimize v. 

2 0 A spectrally constrained impulse shortening filter having an impulse response, g(n), may 

be employed in the receiver to minimize the length of the cyclic prefix by decreasing the EDIR 
of the effective communication channel, which includes the transmit and receive filters, the 
impulse shortening filter, and the physical transmission channel. The use of an impulse 
shortening filter is referred to as time domain equalization. Decreasing the EDIR allows a 

2 5 shorter cyclic prefix to be used without increasing ISI. 

Fig. 2 is a plot of the EDIR for a DMT test configuration having a communication 
channel that is 4500 feet in length and operates at a sampling frequency of 1 1 .04 MHz (test 
loop 4, as described in "Very-high Speed Digital Subscriber Lines: System Requirements," 
Technical Report Tl El .4/97- 1 3 1 R 1 ANSI 1998). The EDIR includes the effects of a transmit 

3 0 filter, the communication channel and a receive filter. Fig. 3 shows the impulse response as 

shortened by the addition of an impulse shortening filter. 

The impulse shortening filter is selected so that a significant part of the energy of the 
joint impulse response of the filter and the effective communication channel, g(n)*h(n), is 
confined to a region that is shorter in length than the length of the cyclic prefix. Some prior 
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algorithms for computing g(n) considered only shortening the EDIR and did not consider the 
spectral characteristics of the resulting impulse shortening filter. Such filters often had deep 
nulls in some frequency bands, which rendered some of the corresponding subchannels useless. 

Since increasing the length of the cyclic prefix reduces system efficiency, the receiver 
5 may dynamically compute an optimal length for the cyclic prefix and may send that information 
to the transmitter. For example, the receiver may compute a set of impulse responses for the 
impulse shortening filter based on a set of predetermined cyclic prefix lengths. The receiver then 
computes the system throughput for each particular cyclic prefix length. The length that 
maximizes system throughput is selected and the result of that selection is communicated to the 
1 0 transmitter. The transmitter then operates using the selected cyclic prefix length. 

To avoid the possible attenuation of frequency bands, the spectral response of the impulse 
shortening filter is further required to have a spectral response that |G(cd)|, meets a specified 
spectral constraint. A spectral constraint of the form |G(© )//(©))> r , where r is a threshold, is 
sufficient to avoid nulls in the frequency response of the impulse shortening filter. However, it is 
15 possible to compute a spectral constraint or desired spectral response, |GX<*>)|, that provides 

additional performance improvements, such as reducing noise bleeding between subchannels. A 
spectrally constrained impulse filter is configured to have a spectral response that approximates 
the desired spectral response. 

As shown in Fig. 4, the spectrally constrained impulse shortening filter (SCISF) 90 may 
2 0 be implemented as a time domain digital filter, which has a digital filter structure 92 with a 

number of taps 94 or filter coefficient inputs for adjusting the filter response. The coefficients 
may be computed and supplied to the taps by a digital signal processor (DSP). Alternatively, the 
SCISF may be implemented entirely in software, i.e., within a DSP. 

A desired spectral response may be applied to a received signal using a filter that is 

2 5 separate from the impulse shortening or time domain equalization (TEQ) filter. For example, an 

analog filter may be placed prior to the A/D converter. However, the adjustability of such a filter 
would be limited. As a further example, a digital filter could be added prior to the TEQ filter. 
Both of these configurations suffer the disadvantage that the TEQ filter may distort the desired 
spectral characteristics of the added filter. A filter also might be positioned after the TEQ filter, 

3 0 which would reduce noise bleeding, but might reduce the impulse shortening provided by the 

TEQ filter. Accordingly, the SCISF integrates the TEQ (i.e., impulse shortening) function with 
the desired spectral response in a single filter. 

In summary, the filter characteristic, g{n), of the SCISF satisfies two conditions. First, 
the effective length of the convolution of the filter characteristic with the impulse response of the 
35 communication channel, g(n)*h(n\ is less than a target length. Second, the cost function (error 
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function) between the desired spectral response, G^o), and the actual filter spectral response, 
G(o>), is minimized. 

The desired spectral response is an ideal filter characteristic that is selected to maximize 
the data bit throughput in the subchannels of a DMT system by reducing the impact of noise. 
5 There are many sources for noise in a DMT communication system, such as near-end cross-talk 
(NEXT), radio frequency interference (RFI) and noise generated in the communication channel 
(white noise). As shown in Figs. 5 and 6, the noise spectral density is generally not uniform 
across the frequency band of the communication system. This non-uniformity contributes to the 
problem of noise bleeding, in which noise in one frequency band interfering with a signal in 
1 0 another frequency band. 

In general, noise bleeding is caused by sidelobes of filters in a DFT. Fig. 7 shows the 
filter response of a DFT for one frequency band or bin (i.e., bin 128). The first sidelobes (96) are 
only 13 dB below the main lobe (98). Therefore, noise located outside of bin 128, but within the 
first sidelobe of bin 128, i.e., approximately mid- way between bins 126 and 127, would appear in 
15 bin 128 with an attenuation of only 13 dB. Consequently, noisy subchannels in a DMT system 
may degrade the performance of non-noisy subchannels. 

The desired spectral response is essentially a spectral constraint that attenuates noisy 
channels more than non-noisy channels. The signal and noise in the noisy channels are reduced 
equally, so the attenuation does not affect the signal-to-noise ratio in these channels. However, 

2 0 because the absolute noise level in the noisy channels is reduced, there is less noise available to 

enter the sidelobes of the non-noisy channels. Hence, the noise bleeding problem is minimized. 

To determine the desired spectral response, the noise power spectral density (noise PSD) 
at the receiver must be known. The noise PSD may be determined, for example, by performing a 
periodigram on received data. This measurement is more complicated if a transmitter is 
25 transmitting, since the transmitted signal must be separated from the noise measurement. The 

noise PSD is determined by: (i) slicing the received constellation of subcarriers after the DFT to 
determine the nearest valid constellation point; (ii) determining an error signal based on the 
difference between the received constellation point and the valid constellation point; (iii) 
performing an IDFT on the error signal; and (iv) generating a periodigram (with windowing) 

3 0 from the error signals. The noise PSD may then be determined from the periodigram. 

An example of a noise PSD characteristic for a DMT communication system is shown in 
Fig. 5 (test loop 4, as described in "Very-high Speed Digital Subscriber Lines: System 
Requirements," Technical Report T1E1 .4/97-1 3 1R 1 , ANSI 1998). The transmit signal power is 
measured at the output of the transmit filter in the transmitter. The signal and noise PSD plots 
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shown in Fig. 5 are measured at the input of the A/D converter in the receiver, which is prior to 
the SCISF. 

Measured noise PSD is used by a digital signal processor (DSP) to compute a desired 
spectral response, G^co), usin 8 * e algorithm described below. Alternatively, the inverse of the 
5 noise PSD may be used as an approximation for the desired spectral response. A spectral 

response, G(<»), is then determined for the SCISF that minimizes the error between the spectral 
response of the SCISF and the desired spectral response. A set of filter coefficients may then be 
generated to configure the SCISF to the determined characteristic. These calculations may be 
done periodically to adjust the performance of the communication system. Frequency domain 
10 equalization coefficients and symbol synchronization may also be adjusted based on these 
calculations. 

Fig. 8 is a plot of the desired spectral response of the SCISF, G d (eo), versus the actual 
frequency response, G(co). The difference between the responses is only a few dB. Fig. 6 shows 
the signal and noise PSD at the output of the SCISF. 

15 Fig. 9 shows the dramatic effect of the SCISF on system performance. Without the 

SCISF (i.e., using a filter that provides only impulse shortening), the signal-to-noise ratio (SNR) 
decreases significantly at the output of the Fourier transform (i.e., FFT or DFT) to less than about 
7 dB. This decrease is due, in large part, to noise bleeding caused by the sidelobes of the Fourier 
transform. By contrast, with the SCISF, the SNR at the output of the Fourier transform tracks 

2 0 the SNR at the output of the SCISF within a few dB. Overall, the SCISF provides an 
improvement in SNR. 

The desired spectral response of the SCISF, G^a)) and the actual frequency response, 
G(co), are derived from an energy constraint for the SCISF, i.e., the desired spectral response 
must localize the energy of the effective impulse response within a desired frequency band. The 

2 5 energy constraint is combined with a desired spectral response based on the measured noise 

power spectral density. The resulting cost function (or error function) is then minimized to 
obtain a practical filter characteristic for the SCISF. This process is presented in detail below. 

In the following derivation, all vectors are column vectors by default. Vectors are denoted 
by bold lower case letters (e.g., t). The size m t of a vector t, is written t(m,) . The components 

3 0 of a vector are denoted by lower case letters, e.g., t(m,) = [r 0 t mf J. The convolution of length 

m t + m h - 1 of vectors t( W/ ) and h(m h ) is denoted t*h. The Toeplitz matrix of a vector 

is written T^ :IJ j , and T^^x = (t * x)^ m: n ^ . The Discrete Time Fourier 
Transform (DTFT) of a vector t may be represented as: 
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7>) = S l„e- Jam . (1) 

The EDIR of the channel and the impulse response of the SCISF are expressed as h(m h ) 
and g(m g ) 9 respectively. The cyclic prefix has a length of m c per symbol. The energy of the 
effective impulse response must be localized to a contiguous region of target length m h where 
mfZm^ while satisfying a spectral constraint. The energy criterion to be satisfied by the SCISF 
may be written as: 

h ^<— ^>L> a (2) 



for some Q < m < m g + m h ~2 and some 0 < a < l . Defining the set S m as 



s = < 



l|2 



g€R ^ : .r8*^)(^-')|| > Q 
life "Ml 



(3) 



1 0 the impulse response of the SCISF must belong to S m for some m. 

Let a>j be the location of the sub-carriers in the frequency domain. The spectral 

constraint may be applied by selecting g e S m , so that the cost function 

4g)= Z^M-g^,)) 2 (4) 

1=1 

is minimized for the desired spectral response G d (co). 
15 Typically, this optimization would have to be repeated for all possible m to select the 

filter impulse response, g, that achieves the lowest possible value of J. However, as discussed 
below, optimization may be limited to a few well chosen values of m. 

The determination of the filter impulse response may be done in two stages. First, the 
desired magnitude frequency response G d (o>) of the impulse shortening filter g over the bins 

2 0 used by the DMT system is obtained. Second, J is optimized over S m for a specific value of m. 

To determine the desired spectral response G d (co) y it is necessary to use expressions for 
the signal to noise ratios observed in the various frequency bins at the output of the DFT in the 
receiver. The DMT system has M tones (subcarriers), N of which (those from A/j through M 2 ) 
are used, and the communication channel has an analog frequency response H c (f) . The analog 

2 5 noise power spectral density observed at the input of the receiver A/D is S^f) . Prior to 

conversion in the A/D converter, the received analog signal may be filtered by an anti-aliasing 
filter with transfer function H a (f). The EDIR in the absence of the impulse shortening filter is 
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h(n). After the A/D converter, the signal is fed into the impulse shortening filter with an impulse 
response of g(n). The impulse shortening filter ensures that (h(n)*g(n)) is shorter in length than 
the cyclic prefix. G(a>) is the discrete time Fourier transform of g(n). Under these conditions, 
the expected signal energy ^k) observed in bin k at the output of the length-2Af receiver DFT is 
given by: 

where C, is a constant, 1/7 the sampling frequency and D k the transmitted power in bin k. The 
noise power ^k) in bin k is: 



(6) 



0)=- 



where C 2 is another constant and * denotes the convolution. Assuming that the noise in the 
bands corresponding to unused tones is sufficiently attenuated by the anti-alias filter, ^k) is 
approximately equal to: 



(7) 



where ^n) is defined as: 



P i?r sin 2 (A/fe-A)) l 



(8) 



M l -"M 2 are the used tones and C 3 another constant. Defining x to be the vector of frequency 
magnitudes to be solved for as: 



X = 



b k X k 



(9) 



the SNR in bin k can be seen to be of the form , in which b k are scalars and a k are vectors. 
To determine the desired spectral response, x is chosen to maximize the bit throughput. 



Approximating the capacity of bin k by Vsnr(a) , the optimal spectral profile is obtained by 
minimizing the cost function F 9 where: 
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W-ZJ 5 ?- do) 



k=M x 

This minimization is performed over: 

* = {xeR":||x||==l, Jt f >0, \<i<N (11) 

and can be accomplished by any one of a variety of constrained optimization strategies, as 
5 discussed in Bertsekas, Nonlinear Programming , Athena Scientific, Belmont, MA, 1995. A 
median filter may be applied to the output of the optimization algorithm to smooth the resulting 
desired spectral response. 

Let A = H ( m 0 * mft+In 2) and B = H^ :ro+m/ _ 0 . The energy constraint in equation (3) can be 

written as: 

Matrix A has full column rank, since it corresponds to a full convolution, so K A ApJ A is 

invertible. Let R B AB T B . Defining q =?. R° 5 h , where R^ 5 is the square root of the positive 
definite matrix R^, the energy constraint set can be written in terms of q as: 

»j,.^ /'ys .}. 03, 

15 The next step is to reduce the dimensionality, i.e., the number of variables to search over 

in the optimization process. For example, in a video digital subscriber line(VDSL) application, 
an impulse shortening filter having a few hundred taps may be required. Searching over a 
variable space so large is difficult and impractical. Instead, the optimization is performed by 
searching over a cleverly chosen lower-dimensional subset of the variable space. This 

2 0 simplification in the optimization process may be done without significant reduction in the 
performance of the communication system. 

The reduction in dimensionality is accomplished by a transformation of variables. Let 

C 4 R^° 5 R B R^° 5 . Let S = UiU r be the singular- value decomposition of C, where E is a diagonal 
matrix whose (i, /)th element is a i and where the a i are arranged in descending order. If a x < a, 
25 there is no feasible solution corresponding to delay m. If a } > a, let 



U = [U!U 2 } £ = 



0 



Z 2 J 



(14) 
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where U, has size (™ g >™d) and E, has size (m d ,m d ) for some m d . These equations define the 
dimension of the range space of the matrix Uj , over which the search is confined. The 
dimension m d may be chosen either to include all <r i greater than some threshold p, fi < a , or m d 
may be a fixed number. Simulations indicate that m d can be less than m g by more than an order 
of magnitude without significantly affecting communication system performance. 

The dimensionality reduction is achieved by a further transformation of variables: 
q = U,v . The energy constraint set now becomes: 



VA 



using the identity ufcu, = Zj ♦ 



^ jv €R^ : *j ={veR^:v^0, a|v|| 2 - v^v <0 (15) 



1 0 Next, the cost function is expressed in terms of v. For a particular v, the corresponding g 

is given by g = R^ 0,5 U,v, which leads to: 

roG»|)l J"l e- Ja> * **** e" 7 ^" 1 ^ 1 1 

I 5 l-Pg; F- : : i : : | (16) 

Lc?W)J e" 72 "" ... e-^-'^j 

Let D^FAR^ 5 U, . Let D* A real (£>) and D 7 wd imag (D) . D* and D 7 are real matrices of size 
(n, m d ). Let d T Rn ,\<n<N and d] n , l <n £ N be the rows of D* and D 7 , respectively. Then: 

where r w ^d^d^, + d 7 „d£ w . These definitions result in: 

^(v)^y(g(v)) = 2]^/v r r w v-G^(^)1 (17) 



The gradient and Hessian of J Q are 



>)-Zl'-#^|r„v (18) 



"j 

The projection P v (y) of any y eR w d on to Kis defined to be 



20 H^-^h-T^lr^-S^r^rJ 0 9) 
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P v (y) = argmm||y-v| 2 (20) 

T6 develop an algorithm to optimize the cost function J 0 over V 9 an expression must be 
derived for P v (y). There is no closed-form expression for this projection operator, however there 
is a very efficient algorithm to compute it. If y e V , P v (y) = y . If not, the projection onto V is the 
5 same as the projection on to its boundary, defined by: 

boundary (V)A^ eR^ : v * 0, a||v|| 2 - V r L,v = 0 (21) 

The latter may be accomplished using LaGrange multipliers as follows. Defining a modified 
cost function S as: 

S(y)A ||y - vlP + ^(allvll 2 - v^v (22) 

10 The equation v 5 (v)= 0 must be solved on the boundary of V. This reduces to the following 
simultaneous equations: 

[(1+^)1-^,^ = 0 (23) 
a||v|| 2 -v 7 'E 1 v = 0 (24) 
Solving for v in terms of ^ from (23) and substituting in (24), the following equation is obtained 



15 for ^: 



y T diag 



\ a 7 a% J y = o (25) 



Equation (25) may be written as a polynomial equation of order 2m d in ^ The polynomial 
equation must be solved for the real and positive roots, which may be done using one of the 
variety of efficient root- finding algorithms in the literature. One of the real, positive roots must 
2 0 be chosen by enumeration so that, when substituted into the expression for v, it leads to the 

smallest ||y - v|| 2 . 

Since the gradient and Hessian of the cost function J Q and the projection operator onto V 
are available, the optimization can be done very efficiently. For example, the penalty method 
may be used, as described in Bertsekas, Nonlinear Programming , Athena Scientific, Belmont, 
25 MA, 1995, or an iterative strategy consisting of gradient descent followed by projection onto V. 
Both techniques have been tested and perform well. 

The SCISF may be configured, by determining the filter coefficients during a training or 
adaptation period. The SCISF filters the output {y k } of the receiver A/D 80. The coefficients 
are selected using an algorithm that minimizes the squared error between a reference sequence 
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{u k } generated by the receiver and the output of the SCISF {u k } . The SCISF may be a finite 
impulse response (FIR) filter or an infinite impulse response (IIR) filter. The SCISF may be 
trained following activation of the communication system or periodically during operation of the 
system to compensate for variations in the channel noise profile. 

The training may be performed using a variation of one of the classical adaptive 
algorithms, such as least mean squares (LMS), normalized LMS, or recursive least squares 
(RLS). For example, the following algorithm is a version of the normalized LMS algorithm in 
which b 0 ,... y bx and a v ... y a N are the FIR and IIR parts of a SCISF having an impulse response, g. 
The z-transform G(z) is: 

G(2) , y^;'— v-; (26) 

1- a^z 1 a N z /v 

The adaptation of coefficients a/ and 6/ is defined in the following equations, in which ai{k) and 
bi(k) are the values of these coefficients during the iteration. The parameter ^ is a 
1 5 predetermined constant with a value of 0.4. 



10 



i=o y=i 

20 a k =[b 0 (k\...,b N (k)] (28) 

A-foWU^M] (29) 

d* =[y*ov.i.....jw] (30) 

c* =["*., ,"t- 2 ,..., u k _ N ] (31) 

(32) 



25 
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(33) 




(34) 



5 



In a first embodiment, the coefficients of the SCISF are determined during an initial 



training period following the activation of the communication system using the LMS algorithm 
described above. A predetermined sequence of bits x k is input to the transmitter 12. The 
sequence of bits results in a sequence of real numbers {x k } at the input of the D/A 60. As shown 
in Figs. 1 and 10, the transmitted signal is filtered and noise-corrupted by the transmission 

1 0 channel 70, resulting in a received sequence {y k } at the output of the A/D 80 in the receiver 14. 
The SCISF 90 filters and transforms the received sequence {y k } into an output sequence {x k } . 
The output sequence {x k } is compared (in a signal comparator 205) to the predetermined 
sequence x k9 which is stored in memory 210 in the receiver. The comparison results in an error 
signal e k that is input to the LMS algorithm processor 215 along with sequence {x k } . 

15 The training process determines coefficients for the SCISF so that the output {x k } 

matches the predetermined sequence {x k } as closely as possible in a least squares sense, i.e., the 
mean square error between the output and the predetermined sequence is minimized. During the 
training process, the coefficients of the SCISF converge to values that enable the SCISF to 
reduce ISI and additive noise. The resulting SCISF matches p^(co) in the frequency domain in a 

2 0 least squares sense, where: 



In the equation above, S x (o>) is the power-spectral density at the input of the transmitter D/A 60, 



(35) 



25 S n (co) is the power spectral density of the additive noise at the output of the A/D 80, and H(co) is 
the frequency response of the effective discrete-time impulse response (EDIR) of the 
transmission channel 70, transmit filter 65, and receive filter 75 measured between the input of 
the transmitter D/A 60 and the output of the receiver A/D 80. 
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Upon completion of the initial training of the SCISF, the SCISF coefficients are fixed and 
the frequency domain equalizer (FEQ 130) of the receiver is trained using standard techniques 
for DMT receivers. Following the training of the FEQ, the SCISF can be periodically trained or 
adapted during operation of the communication system. Since it is not efficient to repeatedly 
5 transmit a predetermined bit sequence during operation of the communication system, the 
periodic training process uses transmitted communication data to generate the reference and 
output sequences, as described below. 

During operation of the communication system, a sequence of communication data bits x k 
is input to the transmitter 12 (see Fig. 1). The sequence of data bits results in a sequence of real 
1 0 numbers {x k } at the input of the D/A 60. As shown in Figs. 1 and 1 1 , the transmitted signal is 
filtered and noise-corrupted by the transmission channel 70, resulting in a received sequence 
{y k } at the output of the A/D 80 in the receiver 14. The SCISF 90 filters and transforms the 
received sequence {y k } into an output sequence {x\ } . 

The received sequence {y k } is also input to a delay 305 and then to a secondary 
15 SCISF 300 that has the same coefficients as the primary SCISF 90 following the initial training. 
Such a configuration allows periodic training to be performed without disruption to the operation 
of the communication system. The secondary SCISF 300 is periodically or continuously trained 
using an algorithm similar to that used for the initial training. The new coefficients of the 
secondary SCISF 300 are periodically copied to the primary SCISF 90. 

2 0 In an initial training process, the output sequence {x k } of the secondary SCISF 300 

would be compared (in a signal comparator 205) to a predetermined sequence x k stored in 
memory. However, as discussed above, a sequence of data communication bits is used as a 
reference sequence for training rather than a predetermined sequence. As such, the receiver must 
have a way of generating a reference signal to compare to the output of the SCISF. 
25 To compute the reference sequence, the receiver essentially duplicates the encoding and 

modulation processes of the transmitter using the output of the decoder 140. Because the initial 
training has already been performed, the SCISF 90 output {x\ } matches the predetermined 
sequence {x k } closely and ISI and additive noise are minimized. Hence, the data output of the 
decoder 140 closely matches the transmitted sequence of communication data bits x k . The data 

3 0 bits are input to an encoder 320, an IDFT 330, a parallel to serial converter 340, and a prefix 

adder 350 similar to those in the transmitter. The output sequence {x k } of this chain is input to 
the LMS algorithm processor 215 and used as a reference sequence in the training algorithm. 
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The output sequence {x k } of the secondary SCISF 300 is compared (in a signal 
comparator 205) to the reference sequence {x k }, which is output by the encoding/modulation 
chain (320, 330, 340 and 350). As noted above, the received sequence {y k } passes through a 
delay 305 before being input to the secondary SCISF 300. The delay 305 compensates for the 
5 processing delay in the demodulation/decoding chain and the encoding/modulation chain. The 
comparison results in an error signal e k that is input to the LMS algorithm processor 215. The 
training process determines coefficients for the secondary SCISF 300 so that the output {x k } 
matches the reference sequence {x k } as closely as possible in a least squares sense, i.e., the mean 
square error between the output and the reference sequence is minimized. Periodically, the 

1 0 coefficients of the secondary SCISF 300 are copied to the primary SCISF 90. 

Alternatively, as shown in Fig. 12, the periodic training may be performed with a single 
SCISF 90. In this configuration, a received sequence {y k } is output by the A/D 80 in the 
receiver 14. The SCISF 90 filters and transforms the received sequence {y k } into an output 
sequence {x\ } . The received sequence {y k } is also input to a delay 305. After the received 

15 sequence {y k } passes through the SCISF 90, a data switch 360 is changed from position A to 

position B, allowing the delayed received sequence to make a second pass through the SCISF 90. 
An output switch 370 also may be opened, so that data is not output during the training process. 
In addition, the SCISF coefficients are controlled by the LMS algorithm during the training 
process. 

2 0 A reference sequence is computed as in the configuration of Fig. 1 1 . The data bits are 

input to an encoder 320, an IDFT 330, a parallel to serial converter 340, and a prefix adder 350. 
The resulting reference sequence {x k } is input to the LMS algorithm processor. 

The output sequence {x k } of the second pass through the SCISF 90 is compared (in 
signal comparator 205) to the reference sequence. As noted above, the received sequence {y k } 

2 5 passes through a delay 305 before being input to the SCISF 90 for the second pass. The delay 

305 compensates for the processing delay in the demodulation/decoding chain and the 
encoding/modulation chain. The comparison results in an error signal e k that is input to the LMS 
algorithm processor 215. The training process determines coefficients for the SCISF 90 so that 
the output {x k } matches the reference sequence {x k } as closely as possible in a least squares 

3 0 sense, i.e., the mean square error between the output and the reference sequence is minimized. 
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The coefficients of the SCISF 90 then are updated to the coefficients determined in the training 
process. 

In a second embodiment, the SCISF 90 coefficients are chosen so that the frequency 
response of the SCISF matches a desired spectral response G d (co) that seeks to minimize the 
5 effects of noise-bleeding and maximize system bit throughput. The desired spectral response 
G d {co) is determined based on the signal-to-noise ratios observed in the various frequency bins 
of the DFT 120 in the receiver. 

For example, an OFDM system may have M tones, N of which (m, through m N ) are used. 
The system operates over a channel with analog frequency response H c (/) . Referring again to 
1 0 Fig. 1 , the analog noise power spectral density at the input of the receiver A/D 80 is S u (f) . 

Prior to receiver A/D 80, the received analog signal may be filtered by an anti-aliasing filter (i.e., 
receive filter 75) having a transfer function H a (J). The effective discrete-time impulse response 
(EDIR) of the transmission channel of the OFDM system (including the transmit filter 65 and 
receive filter 75) is h(n). The output of the A/D 80 is input to a SCISF 90 having an impulse 
1 5 response g(n). G((o) is the spectral response corresponding to g(ri). 

The expected signal energy ^k) observed in frequency bin k at the output of the DFT 120, 
which has a length of NM 9 is: 

^Ai^WSt- ™ 

20 

where Q is a constant, 1/ris the sampling frequency and £fc is the transmitted power in 
frequency bin k. The noise power r|(£) in bin k is: 



25 



77(0= c 2 



sin 



sin 



(37) 



~ M 



where C 2 is a constant and * denotes a convolution of the discrete Fourier transforms. If the 
noise in the bands occupied by unused tones is sufficiently attenuated by the anti-alias filter 
(receive filter 75), r\(k) is approximately: 
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where x(n) is defined as: 



sin 2 M^-4 
sin 2 ^^-^]). 



(39) 



m,. . .m^ are the used tones, and C3 is a constant. A vector of frequency magnitudes g is defined 
as: 



10 



(40) 



Y G 

The SNR in frequency bin k is -V"*- for scalars r k and vectors s k . The scalars r t are defined by: 



15 



(41) 



and «(/), the /th component of * A is defined by: 



(42) 



20 



To determine an expression for g that maximizes system bit throughput, the capacity of 
each frequency bin k is approximated by log(l + SNR k ). Accordingly, the optimal spectral profile 
is determined by minimizing the cost function F 9 where: 
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Ffe)=-glogfl + ^l (43 ) 



Since G k = JCr^w^ / Af )| 2 , the minimization of the cost function is performed over all 
positive values of G k , as: 



g op/ =argnunF(g) (44) 



where: 



10 G={geR" :||g| = l,G,>0a</<#} 



(45) 



A variety of constrained optimization strategies may be used to solve the above equations for g opf . 

Once the optimal impulse response g opt and desired spectral response G^ay) (which may 
be expressed as G d {Ktn k I M) for a system having M tones) have been determined, a training 

1 5 process is used to adapt the SCISF 90 so that its impulse response g matches the desired spectral 
response. As shown in Fig. 13, the training process may be generalized as a feedback system. A 
reference sequence x k is input to the system. This corresponds to inputting a predetermined 
reference bit sequence to a transmitter. The reference sequence passes through a transmission 
channel 410 having frequency response H(f) (including the physical transmission channel and the 

2 0 transmit and receive filters). Additive noise r|* from the transmission channel is represented in 
this general model as an external input 420 to the system. The resulting signal y k is input to a 
filter 430 having a frequency response G(f), e.g., a SCISF. The output of the filter 430 is then 
passed to an adaptation processor 440, which computes an error signal based on the feedback 
loop 450 and adapts the filter accordingly. The adaptation processor may, for example, use the 

2 5 LMS algorithm described above. 

The reference sequence x k is also input to the feedback loop 450, which passes the 
reference sequence x k through a scaling filter 460 with frequency characteristic £>(/). The 
frequency characteristic Q(f) of the scaling filter 460 (which may be expressed as a set of 
frequency domain scaling factors Q k ) is determined so that the SCISF adapts to the desired 

3 0 spectral response. The output of the scaling filter 460 is used a reference for the calculation of 

the error signal in the adaptation processor 440, as described above. 
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Using the general feedback system shown in Fig. 13, a SCISF having an impulse 
response g may be trained to minimize the error \\q * x - x * g * h - tj * gf . The resulting filter 
matches Pi(a>) in the frequency domain in a least-squares sense, where: 

*M- ffiTWv (46) 

S x (*>) is the power-spectral density at the input of the system, H(co) is the frequency response of 

the effective discrete-time impulse response (EDIR) of the transmission channel, S n (co) is the 

power spectral density of the additive noise, and Q(a>) is the spectral response of the scaling filter 
1 0 460 having impulse response q. 

The solution for g opt in the equations above specifies only the magnitude of the spectral 
response of the SCISF. If the SCISF is a FIR filter, a linear phase characteristic may be used. If 
the length of the SCISF is n g , the desired values of G(co) for the frequency bins of interest are: 



G d {*M k I M)^g op Xk)e W { J7tM "^ g l \ (47) 

\ 2M ) 



The values Q k are defined by: 



^ k S n {jnk /M)H* {jnk IM) K } 



The values of Q k may be computed during an initial training period and may be periodically 
updated during operation of the communication system. 

As shown in Figs. 14-16, the general feedback training process may be used to perform 
an initial training of a SCISF followed by periodic training analogous to the process described 
2 5 above with respect to Figs. 10-12. One difference between the techniques is that a scaled 
reference signal (x * q) k is used rather than an unsealed reference x k . 

Referring to Fig. 14, to perform the initial training, a predetermined sequence of bits x k is 
input to the transmitter. The transmitted signal is filtered and noise-corrupted by the 
transmission channel, resulting in a received sequence {y k } at the output of the A/D 80 in the 
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receiver. The SCISF 90 filters and transforms the received sequence {y k } into an output 
sequence The output sequence {x k } is compared (in a signal comparator 205) to a scaled 

reference sequence (x*q) k . 

The scaled reference sequence is computed from a copy of the predetermined sequence x k 
5 that is stored in memory 21 0 in the receiver. As a first step, the predetermined sequence is input 
to a serial to parallel converter 510 and a DFT 515. The resulting frequency domain signal is 
input to a scaling filter 520 which applies the set of frequency domain scaling factors Q k that 
causes the SCISF to adapt to the desired spectral response, as discussed above. The scaled signal 
is input to an inverse discrete Fourier transform 330, a parallel to serial converter 340 and a 

1 0 cyclic prefix adder 350, resulting in a scaled reference sequence {x*q) k . The comparison of the 
output sequence {x k } to the scaled reference sequence (x*q) k results in an error signal e k that is 
input to the LMS algorithm processor 215 along with sequence {x k }> Alternatively, a frequency 
domain reference (e.g., a predetermined bit sequence that has been processed by a serial to 
parallel converter and DFT) may be stored in memory in the receiver, which would eliminate the 

1 5 need for the serial to parallel converter and discrete Fourier transform in the feedback loop. 

Following the initial training, the SCISF is periodically trained during operation of the 
communication system. A sequence of communication data bits x k is input to the transmitter. 
Referring to Fig. 15, the transmitted signal is filtered and noise-corrupted by the transmission 
channel, resulting in a received sequence {y k } at the output of the A/D 80 in the receiver. The 

2 0 SCISF 90 filters and transforms the received sequence {y k } into an output sequence {x\ } . 

The received sequence {y k } is also input to a delay 305 and then to a secondary 
SCISF 300 that has the same coefficients as the primary SCISF 90 following the initial training. 
The secondary SCISF 300 provides output sequence {x k } y which is compared to a reference 
sequence during the periodic training process. Such a configuration allows periodic training to 
25 be performed without disruption to the operation of the communication system. The secondary 
SCISF 300 is periodically or continuously trained using an algorithm similar to that used for the 
initial training. The new coefficients of the secondary SCISF 300 are periodically copied to the 
primary SCISF 90. 

To compute the reference sequence, the data output of the decoder 140 is input to an 

3 0 encoder 320. The resulting frequency domain signal is input to a scaling filter 520 which applies 

the set of frequency domain scaling factors Q k that causes the SCISF to adapt to the desired 
spectral response, as discussed above. The scaled signal is input to an inverse discrete Fourier 
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transform 330, a parallel to serial converter 340 and a cyclic prefix adder 350, resulting in a 
scaled reference sequence (x*q) k . The comparison of the output sequence {x k } to the scaled 
reference sequence (x*q) k results in an error signal e k that is input to the LMS algorithm 
processor 215. The training process determines coefficients for the secondary SCISF 300 so that 
5 the output {x k } matches the scaled reference sequence (x*q) k as closely as possible in a least 
squares sense, i.e., the mean square error between the output and the reference sequence is 
minimized. Periodically, the coefficients of the secondary SCISF 300 are copied to the primary 
SCISF 90. 

Alternatively, as shown in Fig. 1 6, the periodic training may be performed with a single 
10 SCISF 90. In this configuration, a received sequence {y k } is output by the A/D 80 in the 
receiver. The SCISF 90 filters and transforms the received sequence {y k } into an output 
sequence {x\} . The received sequence {y k } is also input to a delay. After the received 
sequence {y k } passes through the SCISF 90, a data switch 360 is changed from position A to 
position B, allowing the delayed received sequence to make a second pass through the SCISF 90. 
15 An output switch 370 also may be opened, so that data is not output during the training process. 
In addition, the SCISF coefficients are controlled by the LMS algorithm during the training 
process. 

A reference sequence is computed as in the configuration of Fig. 15. The data output of 
the decoder 140 is input to an encoder 320. The resulting frequency domain signal is input to a 
2 0 scaling filter 520 which applies the set of frequency domain scaling factors Q k that causes the 

SCISF to adapt to the desired spectral response, as discussed above. The scaled signal is input to 
an inverse discrete Fourier transform 330, a parallel to serial converter 340 and a cyclic prefix 
adder 350, resulting in a scaled reference sequence (x*q) k . The scaled reference sequence is input 
to the LMS algorithm processor. 

2 5 The output sequence {x k } of the second pass through the SCISF 90 is compared (in 

signal comparator 205) to the reference sequence (x*q) k . As noted above, the received sequence 
{y k } passes through a delay 305 before being input to the SCISF 90 for the second pass. The 
delay 305 compensates for the processing delay in the demodulation/decoding chain and the 
encoding/modulation chain. The comparison results in an error signal e k that is input to the LMS 

3 0 algorithm processor 215. The training process determines coefficients for the SCISF 90 so that 

the output {x k } matches the scaled reference sequence (x*q) k as closely as possible in a least 
squares sense, i.e., the mean square error between the output and the reference sequence is 
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minimized. The coefficients of the SCISF 90 then are updated to the coefficients determined in 
the training process. 

In a third embodiment, the system dynamically selects the length of the cyclic prefix (CP) 
to maximize data throughput for a communication channel having a particular noise profile. As 
5 discussed above, a CP is added to each symbol prior to transmission through the communication 
channel to reduce the effects of ISI. However, because the CP constitutes redundant data, 
increasing the length of the CP reduces the efficiency of the communication system. Hence, to 
maximize efficiency, the length of the CP must be as short as the noise characteristics of the 
communication channel permit. 
1 0 For a DMT communication system with M tones, the maximum sample rate W 

(samples/second) for a particular channel depends, in part, on the available bandwidth and 
hardware limitations. The sample rate includes communication data and CP bits. For a CP 
length of n C9 the maximum symbol rate (which includes communication data, but not the CP) is 

15 Before determining the optimal CP length, the SCISF should be initially trained to the 

channel. However, a communication system need not have a SCISF to employ the CP 
optimization algorithm. It is noted that the SCISF coefficients determined during the training 
process do not depend on the CP length. The capacity of the sub-channel may be approximated 
as log(l+ SNRj) bits per second, so the number of bits per symbol is £ log(i +SNRj). For a CP 

2 0 length of « c , the maximum bit rate is expressed as a function of the cyclic prefix as: 



The optimal CP length is determined by computing the maximum bit rate for a set of candidate 

2 5 values of CP length and finding the length that maximizes B a {n c ) . 

The signal to noise ratio SNR t of each subchannel is determined by measuring the 
received signal and noise power and computing the ratio of the two. The noise power for the 
I th bin may be measured by transmitting a data communication sequence and computing the 
average of the squares of the errors measured at the output of the receiver DFT. The total 

3 0 received power (signal and noise) S t for the I th bin may be measured by computing the average of 

the squares of the outputs of the receiver DFT. The signal to noise ratio is determined from the 
expression: S t / r, = 1+ SNR f . Since the signal to noise ratio is determined in the receiver, the 



B a (p c )= 



^X,.log(l + flW,) 
2M+ n c 



(49) 
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computed bit rate B a (n c ) must be transmitted back to the transmitter. The transmitter compares 
the bit rate to the values computed for other candidate CP lengths and selects the CP length n c 
with the highest maximum bit rate B a (n c ) . 

Figs. 17-24 show performance simulations for test systems based on system parameters 
5 and test loops described in VDSL Alliance SDMT VDSL Draft Standard Proposal . Technical 
report, ANSI, 1998; and Verv-high-speed digital subscriber lines: System requirement.^ 
T1E1. 4/97-13 1R1, Technical report, ANSI, 1997. The results are for a VDSL system working 
over test loops 2 and 6 of length 4500 feet in the upstream direction. The system has a sampling 
frequency of 1 1 .04 MHz. Noise is generated by near-end cross talk from an interfering ADSL 

1 0 and an interfering HDSL and white noise at a level of -1 40 dBm. The SCISF used in the 

simulations is length- 15 FIR. A version of the Normalized LMS algorithm is used to train the 
SCISF during an initial training period using a predetermined transmitted sequence. 

Figs. 17-20 show the simulated system performance for a communication system having 
the parameters defined for Test Loop 2, which is 4500 feet in length. Fig. 17 shows channel 

15 frequency response with and without a SCISF. The SCISF provides a much more uniform 

frequency response across the frequency band of interest and significantly improves the signal to 
noise ratio (SNR) in the higher frequency bins. Fig. 1 8 is a plot of the error signal 

(10 log \x k -x k \) during the training process. The error decreases rapidly during the first few 
iterations and is nearly converged after only 20-30 iterations. Fig. 19 is a plot of transmitted 
2 0 power spectral density, received power spectral density and the additive noise power spectral 

density over the used subchannels at the output of the receiver A/D. Fig. 20 is a plot of SNR at 
the input to the receiver A/D, which is the maximum attainable SNR. The plot also shows the 
SNR at the output of the receiver DFT without a SCISF and the SNR at the outputs of the 
receiver DFT using an adapted SCISF. 

2 5 Figs. 21-24 show the simulated system performance for a communication system having 

the parameters defined for Test Loop 6, which is 4500 feet in length. Fig. 21 shows channel 
frequency response with and without a SCISF. Fig. 22 is a plot of the error signal 

(10 log \x k -x k \) during the training process. Fig. 23 is a plot of transmitted power spectral 

density, received power spectral density and the additive noise power spectral density over the 

3 0 used subchannels at the output of the receiver A/D. Fig. 24 is a plot of SNR at the input to the 

receiver A/D. The plot also shows the SNR at the output of the receiver DFT without a SCISF 
and the SNR at the outputs of the receiver DFT using an adapted SCISF. 
Other embodiments are within the scope of the following claims. 

BNSDOCID: <WO 9952250A1_I_> 



WO 99/52250 PCT/US99/07422 

30 

What is claimed is: 

1 . A method for equalizing a channel in a multiple carrier communication system, the 
channel being configured to receive a signal and including a spectrally constrained impulse 
shortening filter, the method comprising: 
5 measuring received noise power spectral density; 

computing a desired spectral response based on the measured noise power; 

selecting a frequency response of the spectrally constrained impulse shortening filter 
based on the desired spectral response; and 

filtering the communication signal with the spectrally constrained impulse shortening 

10 filter. 



2. The method of claim 1, wherein the communication system includes a discrete 
Fourier transform and the noise power spectral density is measured at an output of the discrete 
Fourier transform. 

15 

3. The method of claim 1, wherein the communication system includes a discrete 
cosine transform and the noise power spectral density is measured at an output of the discrete 
cosine transform. 

2 0 4. The method of claim 1, wherein the spectrally constrained impulse shortening 

filter is a time domain digital filter. 

5. A method for selecting an impulse response for a spectrally constrained impulse 
shortening filter in a multiple carrier communication system, the method comprising: 

2 5 measuring received noise power spectral density; 

computing a cost function using the noise power, the cost function being dependent 
on the impulse response; 

reducing the dimensionality of a space over which the cost function is defined; and 
minimizing the cost function. 

30 

6. The method of claim 5, wherein the communication system includes a discrete 
Fourier transform and the noise power spectral density is measured at an output of the discrete 
Fourier transform. 



JNSDOCID: <WO 9952250A1J_> 



WO 99/52250 PCT/US99/07422 

31 

7. The method of claim 5, wherein the cost function is used to compute coefficients 
for the spectrally constrained impulse shortening filter. 

8. A method for equalizing a channel in a multiple carrier communication system, the 
5 channel having an impulse response and being configured to receive a signal having a cyclic 

prefix, the method comprising: 

computing a desired spectral response; 

shortening the impulse response of the channel so that a significant part of an energy 
of the impulse response is confined to a region that is shorter than a target length; and 
1 o filtering the signal based on the desired spectral response. 

9. The method of claim 8, wherein the target length is a length of the cyclic prefix. 

10. The method of claim 8, wherein the desired spectral response is computed using 
1 5 measured noise power density. 

11. The method of claim 10, wherein the communication system includes a discrete 
Fourier transform and the noise power spectral density is measured at an output of the discrete 
Fourier transform. 



20 



12. The method of claim 10, wherein the desired spectral response is the inverse of 
the measured noise power spectral density. 



13. The method of claim 8, wherein the filtering step is performed with a filter 

2 5 having a frequency response selected to match the desired spectral response. 

14. The method of claim 8, wherein the shortening of the impulse response is 
performed by a time domain digital filter. 

3 0 15. The method of claim 8, wherein the filtering is performed by a time domain 

digital filter. 

16. The method of claim 8, wherein the shortening of the impulse response and the 
filtering are performed by a time domain digital filter. 
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17. A spectrally constrained impulse shortening filter for a multiple carrier 
communication system, the system being configured to receive a signal and including a channel 
that has an impulse response, the filter comprising: 

an input connected to receive the signal; 
5 a digital filter structure configured to apply a frequency characteristic to the signal, 

the frequency characteristic being determined by filter coefficients; and 
taps connected to receive the filter coefficients, 

wherein the coefficients are selected to shorten the impulse response of the channel so 
that a significant part of an energy of the impulse response is confined to a region that is shorter 
1 0 than a target length and to apply a frequency characteristic to the signal based on a desired 
spectral response. 

1 8. The filter of claim 17, wherein the target length is a length of the cyclic prefix. 

15 19. The filter of claim 17, wherein the desired spectral response is computed from 

measured noise power density. 

20. The filter of claim 19, wherein the noise power spectral density is measured at an 
output of a discrete Fourier transform. 

20 

2 1 . The filter of claim 1 9, wherein the desired spectral response is the inverse of the 
measured noise power spectral density. 

22. A receiver for receiving a multiple carrier signal from a communication channel 

2 5 having an impulse response, the receiver comprising: 

an analog-to-digital converter connected to receive the signal from the 
communication channel; 

a spectrally constrained impulse shortening filter connected to receive the signal from 
the analog-to-digital converter and configured to shorten the impulse response of the channel so 

3 0 that a significant part of an energy of the impulse response is confined to a region that is shorter 

than a target length and to apply a frequency characteristic to the signal based on a desired 
spectral response; 

a discrete Fourier transform connected to receive the output of the spectrally 
constrained impulse shortening filter; and 
35 a decoder connected to receive outputs of the discrete Fourier transform. 
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23. The receiver of claim 22, wherein the desired spectral response is computed from 
measured noise power density. 

5 24. The receiver of claim 23, wherein the noise power spectral density is measured at 

an output of a discrete Fourier transform. 

25. The receiver of claim 23, wherein the desired spectral response is the inverse of 
the measured noise power spectral density. 

10 

26. The receiver of claim 22, wherein the target length is a length of the cyclic prefix. 

27. A modem comprising: 

an encoder connected to receive digital data and configured to output a constellation 
15 of complex values; 

an inverse discrete Fourier transform connected to receive the constellation from the 

encoder; 

a digital-to-analog converter connected to the inverse discrete Fourier transform and 
configured to output a signal to a communication channel; 
20 an analog-to-digital converter configured to receive the signal from the 

communication channel; 

a spectrally constrained impulse shortening filter configured to shorten an impulse 
response of the channel so that a significant part of an energy of the impulse response is confined 
to a region that is shorter than a target length and filter the signal based on a desired spectral 
2 5 response; 

a discrete Fourier transform connected to the filter; and 

a decoder connected to the discrete Fourier transform and configured to output digital 

data. 

30 28. The modem of claim 27, wherein the desired spectral response is computed from 

measured noise power density. 

29. The modem of claim 28, wherein the noise power spectral density is measured at 
an output of the discrete Fourier transform. 

35 
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30. The modem of claim 28, wherein the desired spectral response is the inverse of 
the measured noise power spectral density. 

3 1 . The modem of claim 27, wherein the target length is a length of the cyclic prefix 

5 

32. Software on a processor readable medium comprising instructions for causing a 
processor in a communication system to perform the following operations: 

measure received noise power spectral density; and 

compute a desired spectral response based on the measured noise power spectral 

1 0 density. 

33. The software of claim 32, further comprising instructions for causing a processor 
in a communication system to compute filter coefficients based on the desired spectral response. 

15 34. Software on a processor readable medium comprising instructions for causing a 

processor in a communication system to perform the following operations: 
measure received noise power spectral density; 

compute a cost function using the noise power, the cost function being dependent on 
an impulse response of a spectrally constrained impulse shortening filter; 

2 0 reduce the dimensionality of a space over which the cost function is defined; and 

minimize the cost function. 

35. The software of claim 34, wherein the cost function is used to compute 
coefficients for the filter. 

25 

36. A method for adapting a primary impulse shortening filter in a multiple carrier 
communication system, the method comprising: 

providing a secondary impulse shortening filter, 

comparing an output signal of the secondary impulse shortening filter to a reference 

3 0 signal to compute an error signal, 

computing coefficients of the secondary impulse shortening filter in an adaptive 
processor based on the error signal, and 

replacing coefficients of the primary impulse shortening filter with coefficients of the 
secondary impulse shortening filter. 
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37. The method of claim 36, further comprising: 

decoding an output signal of the primary impulse shortening filter to form output 

data, and 

5 encoding the output data to form the reference signal. 



38. The method of claim 37, further comprising applying a discrete Fourier transform 
to the output signal of the primary impulse shortening filter prior to decoding the output signal. 

10 39. The method of claim 37, further comprising applying an inverse discrete Fourier 

transform to the encoded output data in forming the reference signal. 

40. The method of claim 36, further comprising: 

receiving a digital signal from an output of an analog to digital converter, 
15 applying the digital signal to the primary impulse shortening filter, 

delaying the digital signal to produce a delayed digital signal, 
applying the delayed digital signal to the secondary impulse shortening filter, and 
applying the delayed digital signal to the adaptive processor. 

2 0 41 . The method of claim 37, further comprising scaling the encoded output data with 

a set of scaling factors in forming the reference signal. 



42. The method of claim 41, wherein the scaling factors are determined by: 

measuring received noise power spectral density, 
2 5 computing a desired spectral response based on the measured noise power, and 

computing the scaling factors so that the coefficients computed in the adaptive 
processor provide the secondary impulse shortening filter with a spectral response that matches 
the desired spectral response. 



30 43. The method of claim 42, further comprising applying a discrete Fourier transform 

to the output signal of the primary impulse shortening filter prior to decoding the output signal 
and measuring the noise power spectral density at an output of the discrete Fourier transform. 
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44. The method of claim 41, further comprising applying an inverse discrete Fourier 
transform to the scaled, encoded output data. 

45. A method for adapting an impulse shortening filter in a multiple carrier 

5 communication system having a spectrally constrained impulse shortening filter, the method 
comprising: 

comparing an output signal of the spectrally constrained impulse shortening filter to a 
reference signal to compute an error signal, and 

computing coefficients of the spectrally constrained impulse shortening filter in an 
1 0 adaptive processor based on the error signal. 

46. The method of claim 45, wherein the reference signal comprises a predetermined 
signal stored in a memory in the communication system. 

15 47. The method of claim 45, further comprising: 

scaling predetermined reference values with a set of scaling factors to form scaled 

values, and 

applying an inverse discrete Fourier transform to the scaled values to form the 
reference signal. 

20 

48. The method of claim 47, wherein the scaling factors are determined by: 
measuring received noise power spectral density, 

computing a desired spectral response based on the measured noise power, and 
computing the scaling factors so that the coefficients computed in the adaptive 

2 5 processor provide the secondary impulse shortening filter with the desired spectral response. 

49. The method of claim 45, further comprising: 

receiving a data signal from an output of an analog to digital converter, 

applying the data signal to the spectrally constrained impulse shortening filter, and 

3 0 applying the data signal to the adaptive processor. 
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50. A method for adapting an impulse shortening filter in a multiple carrier 
communication system having a spectrally constrained impulse shortening filter, the method 
comprising: 

receiving a digital signal, 
5 applying the digital signal to the spectrally constrained impulse shortening filter to 

obtain a first output signal, 

computing a reference signal from the first output signal, 
delaying the digital signal to produce a delayed digital signal, 
applying the delayed digital signal to the spectrally constrained impulse shortening 
1 0 filter to obtain a second output signal, 

comparing the second output signal to the reference signal to compute an error signal, 

and 

computing coefficients of the spectrally constrained impulse shortening filter in an 
adaptive processor based on the error signal. 

15 

51. The method of claim 50, wherein computing the reference signal comprises: 
decoding the first output signal to form output data, and 

encoding the output data to form the reference signal. 

2 0 52. The method of claim 5 1 , wherein computing the reference signal further 

comprises applying a discrete Fourier transform to the first output signal prior to decoding the 
first output signal. 

53. The method of claim 51, wherein computing the reference signal further 

2 5 comprises applying an inverse discrete Fourier transform to the encoded output data. 

54. The method of claim 51, wherein computing the reference signal further 
comprises scaling the encoded output data with a set of scaling factors. 

3 0 55. The method of claim 54, wherein the scaling factors are determined by: 

measuring received noise power spectral density, 

computing a desired spectral response based on the measured noise power, and 
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computing the scaling factors so that the coefficients computed in the adaptive 
processor provide the spectrally constrained impulse shortening filter with a spectral response 
that matches the desired spectral response. 

5 56. The method of claim 55, further comprising applying a discrete Fourier transform 

to the first output signal prior to decoding the first output signal and measuring the noise power 
spectral density at an output of the discrete Fourier transform. 

57. The method of claim 54, further comprising applying an inverse discrete Fourier 
10 transform to the scaled, encoded output data. 

58. A multiple carrier communication system comprising: 

a primary impulse shortening filter connected to receive a digital signal and to accept 
coefficients, 

15 a secondary impulse shortening filter connected to receive the digital signal, output an 

output signal, and pass coefficients to the primary impulse shortening filter, 

a reference signal generator configured to output a reference signal, 
a comparator connected to compare the output signal and the reference signal and 
output a resulting error signal, and 
2 0 an adaptive processor that computes coefficients for the secondary impulse shortening 

filter based on the error signal. 

59. The multiple carrier communication system of claim 58, further comprising: 
a discrete Fourier transform connected to receive an output signal of the primary 

2 5 impulse shortening filter, and 

a decoder connected to receive the transformed output signal from the discrete Fourier 

transform. 

60. The multiple carrier communication system of claim 59, wherein the reference 

3 0 signal generator comprises an encoder connected to receive output data from the decoder. 
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61 . The multiple carrier communication system of claim 60, wherein the reference 
signal generator further comprises a scaling filter connected to scale the output data from the 
encoder using a set of scaling factors. 

5 62. The multiple carrier communication system of claim 61, wherein the scaling 

factors are determined by: 

measuring received noise power spectral density, 

computing a desired spectral response based on the measured noise power, 
computing the scaling factors so that the coefficients computed in the adaptive 
1 0 processor provide the secondary impulse shortening filter with the desired spectral response. 

63. The multiple carrier communication system of claim 61, further comprising an 
inverse discrete Fourier transform connected to receive the scaled output signal from the scaling 
filter. 

15 

64. A multiple carrier communication system comprising: 

a spectrally constrained impulse shortening filter connected to receive a digital signal 
and to accept coefficients, 

a reference signal generator configured to output a reference signal, 
2 0 a comparator connected to compare the digital signal and the reference signal and 

output a resulting error signal, and 

an adaptive processor that computes coefficients for the spectrally constrained 
impulse shortening filter based on the error signal. 

2 5 65. The multiple carrier communication system of claim 64, further comprising a 

memory that stores the reference signal as a predetermined signal. 

66. The multiple carrier communication system of claim 65, further comprising: 
a discrete Fourier transform connected to receive the reference signal from the 

a scaling filter connected to scale the reference signal using a set of scaling factors, 

an inverse discrete Fourier transform connected to receive the scaled reference signal. 
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67. The multiple carrier communication system of claim 66, wherein the scaling 
factors are determined by: 

measuring received noise power spectral density, 
5 computing a desired spectral response based on the measured noise power, 

computing the scaling factors so that the coefficients computed in the adaptive 
processor provide the spectrally constrained impulse shortening filter with the desired spectral 
response. 

10 68. A method for adapting a primary impulse shortening filter in a multiple carrier 

communication system, the method comprising: 

selecting a first cyclic prefix length for the communication system from a set of 
candidate lengths, 

transmitting data having a cyclic prefix of the first selected length to a receiver, 
15 computing a first maximum bit rate value in the receiver for the first selected length, 

and 

sending the first maximum bit rate value to the transmitter. 

69. The method of claim 68, wherein computing the maximum bit rate value 

2 0 comprises measuring received signal power and received noise power of the transmitted data. 

70. The method of claim 68, further comprising: 

selecting a second cyclic prefix length for the communication system from the set of 
candidate lengths, 

2 5 transmitting data having a cyclic prefix of the second selected length to a receiver, 

computing a second maximum bit rate value in the receiver for the second selected 

length, 

sending the second maximum bit rate value to the transmitter, 

comparing the first maximum bit rate value to the second maximum bit rate value, 

3 0 and 

selecting the cyclic prefix length corresponding to the greater of the first and second 
maximum bit rate values. 
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